Josephson effect in CeCoIn^ microbridges as seen via quantum interferometry 
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A superconducting quantum interference device (SQUID) was prepared on a micron-sized single 
crystal using a selected growth domain of a thin film of CeCoIn^ grown by molecular beam epitaxy. 
SQUID voltage oscillations of good quality were obtained as well as interference effects stemming 
from the individual Josephson microbridges. The transport characteristics in the superconducting 
state exhibited several peculiarities which we ascribe to the periodic motion of vortices in the micro- 
bridges. The temperature dependence of the Josephson critical current shows good correspondence 
to the Ambegaokar-Baratoff relation, expected for the ideal Josephson junction. The results indicate 
a promising pathway to identify the type of order parameter in CeCoIn 5 by means of phase-sensitive 
measurements on microbridges. 



I. INTRODUCTION 

Heavy fermion (HF) systems are a class of intermetal- 
lic compounds in which electronic correlation effects are 
particularly pronounced. In HF materials the correla- 
tion effects manifest themselves as a strong increase of 
the effective mass of the charge carriers at low temper- 
atures accompanied with a crossover to the behavior of 
a Landau Fermi liquid^. Hybridization between origi- 
nally localized /-states and the itinerant band states of 
these mainly Ce- and [/-based compounds is causing this 
behavior. In special cases, either by chemical substitu- 
tion or by external stimuli, such as applying pressure or 
a magnetic field, a sub-class of the HF system can be 
tuned into a quantum critical stated In this state, char- 
acterized by the shift of a magnetic phase transition to 
T = OK, magnetic fluctuations exert their influence not 
only at T = OK but also at higher temperatures and 
cause strong deviations from Landau's Fermi liquid be- 
havior. Consequently, these materials are said to behave 
as non-Fermi liquids. In these non-Fermi liquids the low 
lying charge excitations cannot be simply described as 
quasi-particle excitations. Nevertheless, superconducting 
(SC) condensates can form in the non-Fermi liquid state 
which is, e.g., the case for CeCoIn^. The nature of the 
SC state in this HF material is therefore particularly in- 
teresting. Several experiments, such as magnetic field 
dependent heat conductivity^, specific healP and also 
point contact spectroscopy^ have provided evidence for a 
d x 2_ y 2 order parameter. As has been demonstrated for 
selected high temperature superconductors, phase sensi- 
tive experiments employing the Josephson effect can be 
very helpful in elucidating the order parameter type in a 
very direct way. This pathway for HF superconductors 
is, however, difficult to follow. It mainly relies on the 
availability of epitaxial thin films. For selected materials 
efforts in this direction have been successfuP^ and have 
eventually led to well characterized tunnel junctions^HI21. 
For CeCoIn 5 , however, attempts by severals groups have 
so far not led to epitaxial growtH^"2D As a possible 
way out of this difficulty we have previously developed 
a technique for addressing individual growth domains or 
microcrystals in c-axis oriented thin films of CeCoIn§ 



by combining focused ion beam (FIB) etching and fo- 
cused particle beam induced deposition with appropriate 
precursors^HUl. f n this work we take advantage of the 
de- wetting tendency during early stages of CeColn^ thin 
film growth on a-plane sapphire for selecting an individ- 
ual growth domain or microcrystal of about two microm- 
eter lateral dimensions suitable for the fabrication of a 
microbridge-based SQUID. Both, the SQUID loop and 
the individual microbridges exhibit well-defined voltage 
modulations in an external magnetic field. The micro- 
bridges critical dc Josephson current is in good agreement 
with the Ambegaokar-Baratoff (AB) modef^. Further- 
more, in the dynamic resistance we see clear signatures 
of periodic vortex motion in the microbridges. These re- 
sults point towards a very promising approach for phase 
sensitive studies in HF superconductors. 

II. PREPARATIONS 

The CeCoIns thin film was grown on a 10 x 10mm 2 
a— plane a — AI2O3 substrate by molecular beam epi- 
taxy method with the growth details and a typical sur- 
face morphology reported earlier-^. The thin film was 
pre-patterned with conventional ultra violet photolithog- 
raphy and ion beam etching. The micro-structuring (see 
Fig. 1 (a) and (b)) was performed with FIB milling in 
an FEI Nova NanoLab 600 dual beam scanning electron 
microscope (SEM). The outer electrodes (not shown in 
the figure) were metallic ^-composite leads prepared by 
FIB induced deposition using W(CO)q as precursor gas. 
A SEM micrograph of the prepared SQUID is shown in 
Fig. lc and a detailed view of one of the microbridges in 
the inset. The electrical measurements were performed 
in a 3 He cryostat employing a four-probe method and 
a standard differential resistance measurement technique 
using current modulation. 

III. MEASUREMENTS 

The four-probe resistance of the loop was measured 
during cool down at a small constant current /, and 
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FIG. 1. (Color online) (a) FIB milling of the microcrystal 
with its current and voltage leads out of the interconnected 
crystals, (b) FIB milling of the loop center and the bridges, 
(c) SEM micrograph of the such prepared SQUID. The inset 
in (c) shows a detailed view of one of the two bridges. The 
loop dimensions are 2.\[irn x 1.9/im, as a mean value between 
the outer and inner diameters, respectively. The estimated 
microbridges length and width are L m 90nm and width W w 
250nra, respectively. The film thickness is about 300nra. 



its low temperature region is shown in the inset of 
Fig. 2. The complete curve shows the typical be- 
havior for CeCoIns, as was reported for microcrystal 
measurements before^. The residual resistance ratio 
R300K /RaboveT c was found to be 2.9, and is slightly better 
than in the aforementioned work. However, this value is 
still small when compared to the value for bulk crystals, 
which indicates a lower crystal quality for this microcrys- 
tal. The resistivity at 2.5K is also slightly larger than was 
reported for single crystals^. The onset of superconduc- 
tivity is found at 2.0K, somewhat lower than that of 
a bulk crystal. Another smooth transition is found at 
around 1.5K, which we attribute to the Josephson cou- 
pling temperature Tj. This transition we ascribe to the 
SC transition in the bridges, and the local T c being re- 
duced by the FIB milling process. Also, an important 
feature is the finite resistance R Q ss 0.07f2 (also when 
I — > 0) at the lowest temperature, which we explain as 
follows. Thermally activated processes may fundamen- 
tally modify the voltage-to-current V — I relation of a 
Josephson junction. In particular, as was shown by Am- 
begaokar and HalperirPS, there is always a finite resis- 
tance even below the Josephson critical current I c . This 
may be viewed as an overdamped junction, i.e. a junc- 
tion for which the McCumber parameter is small, which 
is expected to be the case for microbridges^Sl. 

The theoretical analysis of bridges whose dimensions 
are large compared to the Ginsburg-Landau (GL) coher- 
ence length £ is difficult. For large ratios of W/£ and 
L/£ the bridge leaves the regime of the 'ideal' Joseph- 
son behavior and enters the Abricosov vortex motion 
regime, as was pointed out by LikhareJ^. As a re- 
sult, essential characteristics of the bridge are modified, 
such as the shape of SQUID voltage modulations or the 
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FIG. 2. Set of measured dynamic resistance curves for the 
SQUID as a function of current. The features on the low- 
est curve are ascribed to the periodic motion of vortices and 
shown with arrows (see text for details). The inset shows the 
measured four-probe resistance during cool down in the low-T 
region. The arrow denotes an additional transition at 1.5K 
associated to the Josephson coupling temperature Tj. Note 
the non-zero resistivity even at the lowest temperatures. 

dependence 7 c (Tp2122l. The dimensions of our bridges 
when compared with £,ceCoi7i 5 ~ 5nnW^ make both ra- 
tios W/S, « 18 and L/\ « 50 large. 

In the main part of Fig. 2 we show a set of measured 
dynamic resistance curves as a function of current plotted 
for selected temperatures. The curves measured above T c 
show an anomalous background resistance, visible up to 
10-fT, out of which evolves the SC state. This nonlinear- 
ity, which we do not attribute to self heating effects, will 
be addressed in a separate publication. Signatures of su- 
perconductivity become evident on top of the background 
anomaly at about 2.0K and become more pronounced as 
the temperature decreases. At temperatures below Tj 
the curves develop additional features, which we marked 
by arrows on the lowest curve. 

The inset in Fig. 3 shows typical measured SQUID 
voltage oscillations as a function of an externally applied 
magnetic field. The oscillations are well reproducible and 
of sine-like shape without asymmetries. The period of 
the measured oscillations corresponds to one flux quan- 
tum $0 = h/2e over the averaged loop area with excellent 
agreement. Occasionally, the oscillations become phase- 
shifted by a random amount. Since the oscillations are a 
measure of the SC phase difference between the two SC 
arms modulated by an external magnetic field, we at- 
tribute these phase slips to randomly trapped flux. The 
main part of Fig. 3 shows a long period modulation super- 
imposed with the SQUID modulations, which appear as 
a band due to the scale. We ascribe them to the Joseph- 
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FIG. 3. Measured voltage oscillations as a function of external 
magnetic field on a large scale. The inset shows the small scale 
SQUID oscillations. Note in the main part that due to the 
scale the SQUID oscillations appear as a band of values on 
top of the slow modulation. 



son modulations of the two individual bridges. These os- 
cillations are also affected by random phase slips. Their 
perturbed periodicity might be caused by interference ef- 
fects between both bridges, since they are coupled by the 
SC wave function in the arms. A rough estimate of the 
period for these oscillations gives lOmT, and according 
to s = &o/B implies that the effective dimensions of each 
bridge extend into each SC bank by about 200nm, which 
compares favorably with the London penetration depth 
X L for CeCoIn^. 



IV. ANALYSIS 

The observed strong manifestation of the Josephson 
effect enables us to analyze the peculiarities of the trans- 
port characteristics at low temperatures in the frame- 
work of models developed for the Josephson effect in 
microbridges^^Hl. The Aslamazov-Larkin (AL) theory^ 
may be used for the analysis of the V ~ I characteris- 
tics of microbridges of a geometry similar to ours when 
vortex dynamics has to be accounted for, i.e. W £ 
and L > £. In the framework of this model the electrical 
current density across the bridge is non-uniform, and is 
described as j(x) — J j^d\J a 2 — x 2 , where J is the total 
current along the bridge, d is the bridge thickness, a is 
equal to W/2, and x is the coordinate along the width 
of the bridge with x = at the bridge center. It may 
be shown that two vortices of opposite helicity are most 
likely to enter the bridge from opposite edges due to the 
locally highest current density there. 

As is shown in the AL model, the total force acting 
on the n — th vortex is F n — —Fl + F s ± F vv . Here, 



Fl is the Lorentz force due to the interaction with cur- 
rent J, F s is the interaction with the edge, and F vv is 
the vortex-vortex interaction term. At low T a vortex 
remains pinned to the edge due to F s as long as the 
counteracting Lorentz force is Fl < F s . Neglecting the 
F vv term when the number of vortices in the bridge is 
small, one may derive the current at which Fl = F s 
as Jc ~ Jo\/(a/2£), where J = c 2 hd/8e\ 2 L . At this 
and larger currents vortex pairs start a periodic viscous 
motion toward the bridge center with successive annihi- 
lation and new creation of a vortex pair at the edges. 
As the J and Fl increase, more than one vortex pair 
may be present simultaneously within the bridge giving 
a discontinuity on the V — I curve. However, due to 
the associated increase of the F vv , any succeeding cur- 
rent span will differ from the previous one. The equation 
of motion of the n — th vortex is given by x n = F n /rj, 
where rj is the viscosity coefficient. Finally, solving this 
equation one finds the period of moving vortices, which is 
equal to the period T in the Josephson frequency relation 
Tvjj = 2-rrh/T = 2eV, where V is the averaged measured 
voltage. Hence, the electrical current J is related with V. 
According to the AL model discontinuities on the V — I 
curves are periodic in the voltage. One may show that 
this periodicity is given by (SI units): 

_ 1 4tt dh 3 1 
° ~ 2V2 Mo a^X 2 eff e 3 ?? 

where X e ff is the effective penetration depth, which 
we take as Xl, since Xl ~ d and the correction is small. 
We attribute the smoothed peaks marked with arrows in 
Fig. 2 each to a new number of vortex pairs moving simul- 
taneously in the bridge. They appear, as expected, just 
above I c . Three corresponding peaks are well resolved, 
which in our case implies that the strong repulsion be- 
tween neighbored vortices of equal helicity limits consid- 
erably their number. This is very probable because the 
average inter-vortex distance a/3 is comparable with the 
vortex radius. Also, large values of ~ 235nr7p2l for 
CeCoIriz, do increase the inter- vortex repulsive coupling. 
The peaks are not equidistant with current, which is also 
expected from the AL model, but the voltage peaks are 
equidistant with a period of about 125/iV with good ac- 
curacy. This periodicity we ascribe to Vo as predicted 
by the AL model. One may obtain the viscosity co- 
efficient rj from Eq. (I). Using Xl = 235nm, the ap- 
proximate thickness in the bridge as d/2 s» I50nm, and 
£ = IfJwrPlwe find i] sa 5.0 x I0" 15 (AVs 2 m- 2 ). The 
agreement with the Bardeen-Stephen theory^ is good, 
which gives ry w 2.6 x 10~ 15 (^4Vs 2 ra -2 ) using £ = lOnm 
and the normal state resistivity p n = 4.0 x 10 _7 r2m. 

A temperature dependence of the reduced Josephson 
critical current I c for the two bridges is shown in Fig. 4. 
The dots represent measured data and the solid lines are 
theoretical curves for several cases using the measured Tj 
and the normal state resistance R n . It is clear that the 
AB modef^l fits the data well, although R n = has 
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FIG. 4. The reduced Josephson critical current for an indi- 
vidual microbridge plotted versus temperature. Circles rep- 
resent the measured data. Solid curves are calculated for: 
(AB)-Ambegaokar-Baratoff, Kulik-Omery anc huk models for 
diffusive regime (KO-I), and for ballistic regime (KO-II). Tj 
denote the Josephson coupling temperature. The measured 
value of the normal state resistance R n ~ 42Q, while effective 
one used for the AB fit is 68£1. The measured I c at the lowest 
temperature is ~ 5.2/iA 



to be used instead of the measured value of 42f2. This is 
not unexpected. On the one hand, the AB fit tends to 
overestimate the value of I c for L > £, since R n growth 
linearly with L but the GL order parameter in the mid- 
dle of the bridge falls off exponentially. On the other 
hand, the AB model is derived for tunnel junctions, that 
is L <C £ and with no scattering within the barrier. In 
contradistinction, models derived for microconstrictions, 
such as that by Ku lik and Omel'yanchuk, in the dirty 
and clean limfP£l2ZI or Ishii's modeP^for long bridges do 
not fit the data. Additionally, the shape of the voltage 
modulations are not asymmetric, as expected in these 



models. We have to state that the observed I C (T) depen- 
dence does not follow the probably expected behavior 
theoretically predicted for microbridges in the dirty or 
clean limit. These models do not take a possible influ- 
ence of Andreev bound states^ for a d-wave order pa- 
rameter into account, which may be one reason for this 
discrepancy. 

V. CONCLUSIONS 

The clean SQUID characteristics observed in this work 
on selected CeCoIn^ microcrystals allow for two main 
conclusions. First, the fact that the micro-patterning 
approach used here is not detrimental to the supercon- 
ducting state in C'eCoIn^ indicates a certain robustness 
of this state. Considering, as a rule, the high sensitivity 
of the superconducting properties in other heavy-fermion 
superconductors, this may point toward an enhanced sta- 
bility of the superconducting state when formed from a 
non-Fermi liquid. Second, the approach followed here 
provides junctions with rather close-to-ideal characteris- 
tics, which is promising with regard to utilizing other mi- 
crobridge geometries. These geometries, such as recently 
proposed by Gumann and Schopohr-- or also a 7r-junction 
structure, which can e.g. be formed by the combination 
of a conventional SC lead prepared by focused ion beam 
induced deposition using the W(CO)q precursor with a 
CeCoIn§ microcrystal, can provide phase-sensitive in- 
formation about the superconducting state in CeCoIn^. 
Work along these lines is in progress. 
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